Slow failure of quasi-brittle solids 
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A new mesoscopic non-equilibrium thermodynamic approach is developed. The approach is based 
on the thermodynamic identity associated the first and second law of thermodynamics. In the 
framework of the approach different internal dissipative channels of energy are taken in account in 
an explicit form, namely, the thermal channel and channels of defect subsystems. The identity has 
a perfect differential form what permits to introduce an extended non-equilibrium state and use 
the good developed mathematical formalism of equilibrium and non-equilibrium thermodynamics. 
The evolution of non-equilibrium variables of a physical system are described by a Landau-based 
equation set expressed through internal or different kinds of free energy connected by means of the 
Legendre transforms. The accordance between the different kinds of energy is possible owing to 
introduction of some trends into the equation subset described the defect subsystems and having 
a nature of structural viscosity. The possibilities of the approach are illustrated on the example of 
quasibrittle solid damage and failure. Taking into account only one type of defects (viz., microcracks) 
and mechanical parameters in an expansion of free energy down to third powers in relative to 
average energy per microcrack, the description of destruction of quasi-brittle solids during long- 
term loading is considered. The consideration allows to find equilibrium and non-equilibrium values 
of the free energy. A qualitative behavior of the system on parameters of the theory is analyzed. The 
destruction of material is described from the uniform positions, both at uniform tension and uniaxial 
compression. Origins of the high stability of mine workings at small depths and their instability at 
large depths are explained. 
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I. INTRODUCTION 



A problem of solid fracture is one of the most diffi- 
cult challenges for modern science. The fracture is a 
compound multilevel process that depends on many fac- 
tors, such as physical properties of material, its struc- 
tural organization and environment and so on. Variety 
of materials and environments predetermines many frac- 
ture mechanisms that are distributed in the range from 
brittle to plastic failure. These are sudden dynamic de- 
struction during quick loading and a slow accumulation 
of microdefects with their subsequent coalescence during 
a long-term static loading. 

Rocks are one of types of poorly studied solids of pe- 
culiar structural organization. Rock properties to resist 
to long-term destructive effects of rock pressure at large 
depths are important component of stabilization of un- 
derground structures. At the same time rocks demon- 
strate a number of common for all solids features of be- 
havior that can be of certain interest for investigation of 
processes that are initiated by a strong external effect 
on a matter and also for investigation of the resulting 
non-equilibrium states. 
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The governing relaxation mechanism in the case of 
quasi-brittle solids is accumulation and growth of micro- 
cracks, its coalescence and branching and finally destruc- 
tion of a solid [il] , U • Previously this mechanism for the 
case of large deformation was considered of general ther- 
modynamic positions Q , [J] , 0] , 0] , H ■ Applicability of 
equilibrium thermodynamics was limited by intermediate 
time scale less than that of necessary for diffusive heal- 
ing of microcracks. It was expected that a real solid has 
already natural weak points, non-continuities and so on, 
which became more intense due to external force action. 

Rock was considered usually as passive mechanical 
continuum, properties of which are evolved just under the 
influence of active external stresses. Yet, it is clear from 
general considerations that in the process of metamor- 
phism at large depths the rocks reach state what is equi- 
librium in thermodynamic meaning for the given depths. 
In these conditions, the formation of an artificial cavity 
results not only to mechanical disequilibrium, but, what 
is much more essential, to thermodynamic disequilibrium 
of rocks. Actually rocks turn out in a non-equilibrium 
state and relax by means of microcrack accumulation, 
subsequent coalescence of them and macroscopic destruc- 
tion. Depending on the level of accumulated internal 
energy, the autowave processes of different kinds, which 
are typical for an active media, are possible [S;] , including 
generation of circular zonal destruction structures around 
mine workings [9|, [lO|. 

In this work on the basis of the principles of thermo- 



dynamics an attempt to formulate a consistent theory 
is made that takes into account both mentioned factors 
energy accumulated in the form of internal stresses and 
work of external stresses. Areas of increased internal 
stresses make additional contribution to internal energy; 
generation of microcracks also results in growth of inter- 
nal energy of a solid of the external work which is equal 
to a total energy of break bonds. 



II. DEFINING RELATIONS 

For infinitely small changes of external parameters of 
a defect-free homogeneous system during its interaction 
with external bodies a change of internal energy (first law 
of thermodynamics) is equal to 



du = TSs + aijde^j 



TijdSeij, 



(1) 



where u is density of internal energy, J /rrt'; T is temper- 
ature, K] s is entropy density (here from external sources 
only), J/K * m?] aij is the stress tensor, H/w? and ef^, 
e" are the reversible and non- reversible parts of defor- 
mation tensor. Here increments of variables, which are 
state functions, are perfect differentials. In this form the 
law expresses external forced and thermal actions on the 
studied object and it fits for description of both equilib- 
rium and non-equilibrium processes. 

Irreversible part of the external work goes to internal 
heat at the cost of entropy production by internal sources 
and to generation and modification of defect subsystem, 
what may be written as 



(j,j5e\) =T6s" 



xSs' 



N 

y^ fkShk > 0. 
fc=l 



(2) 



Here s" is that part of entropy produced by internal 
sources, which has time to go to equilibrium state dur- 
ing the external act, s' is that part of entropy pro- 
duced by internal sources, which remains yet into the 
non-equilibrium state and it lies ahead to relax in the 
next time, x is the non-equilibrium temperature, s(pk, hk 
are average fc-type defect energy (J) and defect density 
(m~^), respectively, N is number of defect subsystems. 
The relation ([2]) is one of equivalent form of the second 
law of thermodynamics. Substituting it to ([1]) and as- 
sociating equilibrium entropy produced by external and 
internal sources one may deduce next identity 



N 



du = Tds + aij de'^j 



Xds' + y^^ipkdhk. 



(3) 



fc=i 



The relation ^ is combination of the first and second 
law of thermodynamics in the form of identity, but not 
in the known inequality form. It turns out possible to 
write it in this form as all internal channels of energy 
dissipation are known for the given model of a solid. The 
sound property of this identity is that all increments of 



the thermodynamic variables are perfect differentials in 
one or another sense. The first two terms are treated 
in the sense of equilibrium thermodynamics; the second 
two terms describe non-equilibrium processes passing in 
a solid. The temperature and stress tensor may be calcu- 
lated by means of differentiation of internal energy with 
respect to full equilibrium entropy and elastic strain ten- 
sor as it does in the classical thermodynamics. The non- 
equilibrium temperature and defect energy are too calcu- 
lated by means of differentiation of internal energy with 
respect to non-equilibrium entropy and defect densities 



du 
9?' 



fk 



du 
dhk 



(4) 



but in the sense of the extended thermodynamic forces 
in the Landau-based equations. Owing to this, the ex- 
pression ([3]) may be regarded as complete differential 
with respect to variables s, e? , s', and h^. and this vari- 
ables determines some generalized non-equilibrium state 
of the system. What does it give? Now one can use 
the simple mathematics of equilibrium thermodynamics 
for the non-equilibrium case. Time evolution of non- 
equilibrium thermodynamic parameters s', and hk can 
be calculated by means of set of equations each of it is 
Landau-Chalatnikov or Ginzburg-Landau equation type 

0,[ni: 



'dt 



du 



= ls 



ds 



7 = 7s'X, 



dhk , du _ _ 

-Qf = ^hA-QJ^ - Vk) = IhA'Pk - <Pk)- 



(5) 



These equations differ from the classical prototype. 
First, here the internal energy is used as a more fun- 
damental quantity entered to the first law of thermody- 
namics. Second, in the right parts of equations the sign 
plus stands it means that the stable stationary states are 
placed near a maximum of the internal energy. This pe- 
culiarity is conditioned by energetic pumping in the cost 
of irreversible work. The third difference is presence of 
viscosity shift Ipj. in the equation describing defect sub- 
system. In consequence of the shift the stable stationary 
states do not coincide exactly with the maximums of the 
internal energy, but coincide with some points in which 
the curve of the internal energy has finite positive slope 
of tangents. This shift has same physical sense as order 
parameter in the case of a supercooled liquid and plays 
an important role for construction of a consistent set of 
non-equilibrium thermodynamic potentials connected by 
means of the Lcgendre transformations. 

So as number of conjugate pairs of thermodynamic 
variables increases in comparison with equilibrium case 
the number of the Legendre transformations increases 
too. For example, if one carries out the Legendre trans- 
formation f — u — Ts the classic free energy is resulted. 
Difference of it from classic one is that this type of free 
energy has non-equilibrium nature. The stable station- 
ary states of a system in reference to defect subsystem 



are placed near the maximums of the free energy same 
as for internal energy. 

If one carries out the Legendre transformation in the 
form 



N 



Ts 



y^^i-pkhk 



(6) 



fc=i 



the complete free energy is resulted. Despite to the above 
the stable stationary states of the system in reference to 
defect subsystem are placed near the minimums of this 
free energy. Its total differential can be written as 



df 



-sdT 



O'ijQ'Oij 



Xds' 



N 

^hkdipk- 



(7) 



From this expression one notices that the complete free 
energy is function of such own arguments as T, e? , s' and 
ifk- Time evolution of non-equilibrium thermodynamic 
parameters s' and (pk can be calculated by means of set 
of equations which is similar to ([5]) 



(8) 



The non-equilibrium temperature and defect density 
are calculated by means of differentiation of the complete 
free energy with respect to non-equilibrium entropy and 
defect energy 





ds' df 

dt ~^'' ds' ~^''^' 


dt 


^IvAq hk)^lvk{hk hk). 



X 



5/ 
ds'' 



dipk 



(9) 



These relations and relations (H]) can be interpreted 
as the state equations. Ignoring the thermal effects and 
supposing that the drive parameter ef is a constant the 
state equation can be written in two equivalent forms: 



<^fc 



du 
dhk' 



dipk 



(10) 



III. POWER SERIES ANALYSIS 



Viewing only one type of defect power series expan- 
sions of the internal and free energies in terms of its own 
arguments can be written in the above approximation as 



u = uo + Voh - -Vih'^ + ■r'^2/1^ 



1 



ip,h\ (11) 



/ = /o - hotp + -^hup^ - -h2ip^ + -rhs^'^- (12) 

For simplicity sake, only one type of defects is consid- 
ered here. The coefficients ipo, ipi, ip2 and ips characterize 



energetic processes at microscopic levels, the coefBcients 
/iQ, hi, /i2 and hs characterize dependence of defect den- 
sity from processes at microscopic levels. Different pow- 
ers in the expansions describe thermodynamic processes 
of different levels. The main process is described by the 
lowest power; the processes corrected to it are described 
by higher powers of the expansions. Alternation of the 
signs is due to the Le Chatelier principle the processes 
of each next level point in opposite direction to the pro- 
cesses of the previous level. The state equations for each 
case have the form: 



•9" , ,2 , s 

1^ = — - = (po - Vih + tp2h - tpsh , 
dh 



h 



dip 



Hq - hup + /i2(^"' - h-^if' 



(13) 



(14) 



Let us notice that in the linear (quadratic for energy) 
approximation the dependence between thermodynamic 
variables ip and h is biunique in the both cases. Taking 
in account that in this approximation the both equations 
must lead to the same physical result the connection be- 
tween expansion coefficients is next: 



ho 
hi 



1 
hi 



ho = 



'£o_ 
'Pi 



hi = 



Pi 



(15) 



(16) 



Calculation of the internal and free energy and the 
state equation at model parameters of uq = 0.114 • 
106 J/m3, Lpo ^ 0.113 • 10-lJ/m^ <pi = 0.0.976 ■ 
10~^^J/to is presented in the fig. [T^. The dimension 
of the variable h is m~^. If in accordance with the state 
equation (|13p . (|14p written in the linear approximation 
one makes transition from variable h to variable ip. then 
sequence of extremes are changed, the maximum of the 
internal energy (at zero value of ip) is placed in the left 
side, the minimum of free energy is placed in the right 
side. 

From the fig. [T^ one notices that the maximum of the 
internal energy does not coincide with the minimum of 
the free energy. The standard procedure of minimization 
of free energy has led to the zero value of the variable 
h as in the same time the analogous procedure of max- 
imization of internal energy has led to the zero value 
of the variable ip suggesting that there is some conflict 
between descriptions in terms of the internal and free en- 
ergy. It is possible to eliminate this conflict offsetting the 
extremums of potentials (energies) to meet one another 
in such away that they both turn out in the same point 
satisfying to the state equation. In the terms of the ini- 
tial potentials (energies) it means the stable point is that 
in which the slopes of tangents to the internal and free 
energy curves are not equal to zero and together they are 
satisfied to the state equation. In fact, an arbitrary point 



■6 

s 
o 

s" 
o 
■a 

B 
3 

(D 



So 

43 



I \ 


1 J 


1 


\i \ 


2/ 




\ 
\ 




i 
i 
i 
i 

i 
i 


i\ 




i/ 


\^ 


\. 


/ 





■6 

s 
o 

s" 
o 
■a 

B 
3 

(D 



43 



B 

s 



1 


1 


f 


\ 




/l 






i 
i 
i 
i 






i 
i 


l/ 


i 
^^^- + - 

1 


\ 



a 



0.2 

Defect density, cond. units 



0.2 

Defect density, cond. units 



FIG. 1: Dependence of the thermodynamic energies versus defect density: a ■ 
2 free energy; 3 state equation (p = tp{h) 



initial ones; b - shifted ones; 1 internal energy; 



between zero and h^ax satisfies to this condition, conse- 
quently the position of the point must be determined 
from experimental data. The plots of shifted thermody- 
namic potentials for value h — O.lm^^ are presented in 
the fig. \Tjp. The value of Ip determined from the state 
equation is 0.15 • Jm~^. 



From the figure we notice that the curve of internal 
energy is shifted a little to the left side, and the curve 
of free energy is essentially shifted in the right side so as 
the extremums of both energies are coincided. Thus, a 
system during its evolution tends to common stable state 
independently of what kind of energy is used for descrip- 
tion. In consequently of it the evolution of a system can 
be described by two kinds of Landau-Khalatnikov equa- 
tions expressed in terms of internal energy ([5]) or free one 



All along, it is possible to determine the kinetic co- 
efficients 7/i and 7,^ so that the equations ([5]) and ([5]) 
at a arbitrary time will describe one and the same state 
associated with thermodynamic pair h and if connected 
by the state equation P^ or (|14p. The evolution is fin- 
ished with achievement of one and the same system state 
associated with pair h and Tp. 



IV. FAILURE OF QUASIBRITTLE MATERIALS 

A. General 

In a special case of isothermal processes, which is just 
considered here, valid is 



df 



T.jdel^ 



N 

E 

k=l 



hkdifk, 



(17) 



so / = fi^ij, Vk) If explicit expression of the free energy 
density dependence of its arguments is known, then vari- 
ables Uij and hk can be determined by simple differentia- 
tion. As determination of the explicit dependence for the 
density of the free energy is intractable problem, we use 
a standard method for such cases, expansion of the free 
energy into a series up to cubic power of its arguments: 



/(4,'^) = /o 



hoV+ -^hi(f^ 



"^u 3 



1 



+2^4)' + M4) 



+ Xp{e%f +-pp{sl^f 



gv^t 



eV'^e?^. 



(18) 



Here we consider only one type of defect, which is main 
for quasibrittle failure, namely, micro-cracks. Then h 
is micro-crack density and Lp is average energy per one 
crack. During an external mechanical action the micro- 
cracks are produced, collected and merged into macro- 
scopic clusters such as main-cracks and so on. 



Corresponding state equations have the form: 

h = ho — hiip + /i2(^^ + .gSjj + 

Gi-j = XsiiSij + 2^eij - gipSij (19) 

-XipSiiSij - 2jLipeij + eip^5ij. 

The evolution equation in an exphcit form now is 

-^ = Xv{ho -h- hi<fi + /i2<y3^ + 9£u + 

+ lx{£u)^+neW'-^e^£H)- (20) 

The positive terms represent microscopic processes as 
a contributory factor for increasing of average energy of 
separated micro-cracks; the negative terms are same for 
decreasing energy. The constants ho and h are constant 
sink and source for average micro-crack energy. Micro- 
scopic sense of the constant Hq is foUowed; it represents 
that part of defectiveness of a deeper structural level 
which provokes activation of new micro-cracks. Follow- 
ing from it increasing of micro-crack density reduces the 
elastic property of the material and, as consequence, the 
average energy of a micro-crack. Constant h represents 
that part of defectiveness of the deeper structural level 
which provokes micro-crack healing, what leads to de- 
creasing of defect density and increasing of the average 
energy of the remaining micro-cracks. Regeneration of 
defects of the deeper structural level is supported by ir- 
reversible work. It is notice that in a particular problem 
the constants ho and h may be integrated in one effective 
constant ho — ho — h. 

The term hitp also represents micro-crack healing but 
connected with surplus energy relaxation of the self 
micro-cracks (self-action). The term h2<p'^ represents 
micro-crack merge. At small ip this process is weak and it 
is easily compensated by the above described processes. 
At large ip micro-crack density increases sharply, in the 
same time the spatial distribution of the micro-cracks is 
not totally independent, but forms chains of macroscopic 
length from merged micro-cracks, that is, new objects of 
a higher structural level. Although energy of the new 
object (macro-crack) is high the part of energy per one 
micro-crack included in it becomes less, what is repre- 
sented by the minus sign before the term. 

In the given approximation free energy curve is a cu- 
bic parabola (Fig. ^. In accordance with a sign of the 
determinant 

D = {h,+ 2eelf) + ih^iho - ^A(4)' - M^)' ) (21) 

three kinds of cubic parabola are possible: with two ex- 
tremes at D > (curve 2) , monotonic with point of in- 
flexion at D = (curve 3), and monotonic without point 
of inflexion at Z? < (curve 4). Which one of the three 
kinds of cubic parabola will be realized in a particular 



case depends on the relation of material parameters and 
the magnitude and character of a stressed state (value of 
tensor efj) 

The equation of state calculated according to (fT9|) is 
represented by the curve 1 in Fig. [21 It is seen that 
at such consideration, the extremes of free energy corre- 
spond to zero values of defect density, i.e. from the point 
of view of mechanics, to a continuous defect-free solid. 
Defect density between extremes has a negative value, 
and it seemingly has no physical meaning. Because of 
this there is something like a gap between the right and 
left branches of the free energy, which prohibits a spon- 
taneous transition from one branch to another. The plot 
of the state equation for the free energy described by the 
curve 3 just touches the abscissa axis and for that de- 
scribed by the curve 4 it is above the abscissa axis. In 
both cases the defect density decreases at first and then 
increases with increase in the average defect energy. It 
implies that due to internal processes the material in the 
given stressed state (at predetermined efj) first tends to 
be consolidated and only then loosened and collapsed. 
Presence of the free energy gap results in that the transi- 
tion from the consolidation stage to the destruction one 
cannot be realized exclusively through the internal pro- 
cesses. In this case, material is stable and will not be 
destroyed during unlimited time. 

The left extreme of the free energy is the minimum 
and, corresponding to it, state is stable. The right ex- 
treme is the maximum and, corresponding to it, state 
is unstable and provokes destruction of a material. The 
free energy difference of these two states represents an 
energy barrier. Existence of such barrier is a prerequisite 
for sustained stability of a material. At large stress the 
state is described by the free energy in the forms 3, 4. 
In these cases, the energy barrier vanishes, and material 
will be in a state of long-time creep-like destruction that 
will finally result in its microscopic destruction. 



B. Numerical calculations 

Let us assume for numerical estimates the following 
set of reference model parameters: /o = 2 • 10^ Jm"^, 
ho = I0^m~^, hi ^ 1.1 • lO^V^im-^^ /i2 = 0.35 ■ 
IQ-^'^ J-^m-^, X = n = 0.208-10"Pa, g ^ 1.5-10"m-3^ 

-3 



A = 0, 71 = 3.5 • lO^'^m--^, e = 5.6 • lO^V-^m-^ In 
all calculations, the main components of the deforma- 
tion tensor are assumed to be constant and equal to 
^11 = ^22 = ^33 — 0.001 that corresponds to uniform 
tension. 

Changes of a form and height of potential barrier with 
variations of the parameters are illustrated in Fig. [3l 
The decrease in the parameter ho in comparison with 
the reference value increases the height of the potential 
barrier (curve 1, Fig. [3^) that is equivalent to increase 
in material strength. This does not contradict common 
sense and satisfies our intuition, namely, the smaller is 
the damage of a material on a deeper structural level, and 
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FIG. 2: Schematic presentation of density of defects and free energy versus average energy of defects: onset of coordinate 
system; 1 equation of state h — h{ip); charts of free energy: 2 with two extremes; 3 with point of inflexion; 4 declining curve 
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FIG. 3: Dependence of free energy versus average energy per microcrack with variations of parameters: a - ho, h - hi, c - h2. 
Numbers mean the following: 2 reference parameters; 1 decrease ones; 3 increase ones 



the more stable it is. The increase in the parameter Hq 
changes the height of the potential barrier sharper than 
in the decrease case (curve 3, Fig. [5^). It is necessary to 
note, that the free energy minimum at large damage on 
a deeper structural level (curve 3, Fig. [3^) has shifted 
noticeably towards large values of (p. It means that in a 
microscopically more damaged material the states with 
large average energy per a microcrack are much more 
preferential from the point of view of thermodynamics. 
Inasmuch as defects of larger sizes have larger energies, 



the instability of material can be connected with achieve- 
ment by microcracks of some critical sizes or critical av- 
erage energy. 

Smaller values of the parameter hi (curve 1, Fig. Wp), 
on the contrary, result in decrease in the potential barrier 
and describe materials with less strength, and its larger 
values represent materials with higher strength (curve 3, 
Fig. [3}d) . This can be explained by the fact that decrease 
in efhcient parameters of material with increase in the 
parameter hi that taking into account interaction of mi- 



crocracks leads to the growth of adaptation capabihties of 
a material, the reduction of internal stress concentration 
that results in growth of its stability. 

Influence of the parameter /12 on stability is outwardly 
the same as that of the parameter hg (compare Fig. [3]and 
Fig. [3J;) . Namely, if it decreases the stability raises (curve 
1, Fig. [3J;) and if it growths the stability reduces (curve 
3, Fig. [3J;). However, this influence is realized according 
to the other mechanism in comparison with case for the 
parameter hg. Namely, with growth /12 the processes of 
coalescence of old microcracks and generation of new ones 
are stimulated to a greater extent that does result in a 
reduction of material strength. 

The dependence on the other parameters of the free 
energy expansion is determined by what power of the 
parameter ip they stand at. All items having the same 
power with respect to tp, , can be united through the 
effective constants of the theory: 



hn 



9£i 



iA(4 



hi + 2e4 



(22) 



The second term in the expression /Iq is linear with 
respect to the reversible deformation. At the tension 
e^j > stability is reduced with the growth of the re- 
versible deformation (Fig. [3^). At the compression, vice 
versa, e^j < and stability is increased with the growth 
of the reversible deformation that is confirmed by every- 
day experience. The last terms in the expression h^ are 
quadratic with respect to deformation. It means that 
at positive parameters A and JI stability of a material is 
equally reduced, both as in case of the tension and as in 
case of the compression. Further, we will assume that 
A << 71, i.e. shear deformations play a major role in the 
destruction of a material. 

In the above examples, the deformation remained a 
constant value. It is necessary to clarify how a system 
described by the theory behaves in a thermodynamic cy- 
cle. Fig. Hlgives results of calculations of uniform tension 
and three-axial unequal-component compression. In the 
first case, the main components of elastic deformation 
tensor are equal to efi = £^2 = £33 — ^ where e is the 
value given in Fig. [4l In the latter case, they are equal 
to ell = 0; £22 — ^7 £33 = 2e. Such conditions are 
specially chosen to describe a situation in which rocks 
are in the vicinity of a free surface of a mine working 
as in a real situation. It is seen that in both cases the 
height of potential barrier decreases during increase in 
load, and at a definite deformation it becomes equal to 
zero. Material loses stability both in the first and in the 
second case. However, instability mechanisms are differ- 
ent. In the first case it is due with long-term ultimate 
tensile strength, in the second one it is due with shear 
strength. However, critical values for a quasibrittle ma- 
terial in the given theory are not introduced explicitly 
as the fundamental strength constants, they are derived 
from the theory. 



V. EXPLANATION OF EXPERIMENTAL 

FACTS AND APPLIED ASPECTS OF THE 

PROBLEM 



The theory presented here can be applied to qualitative 
explanation of phenomenon of rock failure around mine 
workings. It is well known that at small depths mine 
workings have a rather high degree of stability. Around 
identical workings at large depths one can observe phe- 
nomena of anomalously large displacements of rocks in- 
side the free working space that give large reduction in 
the mine working section area and finally lead to the com- 
plete collapse of it Q . According to the theory, at small 
values of deformation tensor e?-, i.e. at small depths (es- 
timates are given below) for rocks surrounding a mine 
working the condition of existence of equilibrium state is 
fulfilled (the case of Z3 < 0) which is described by free 
energy with one minimum (curve 2 in Fig. U). At large 
values of shear deformations (in elastic approach the vol- 
umetric deformation is equal to zero), i.e. at large depths, 
kind of the free energy dependence is changed, and it is 
already described by the curve without extremes (curve 
3 in Fig. [2^). In this case, a system has no equilibrium 
states (at D > 0) and therefore is permanently collapsed. 

The phenomenon described above suggests an elegant 
method to determine an additional estimate coupling 
among the parameters of the theory. Let us assume that 
mine depth H, , starting from which said phenomenon 
begins to manifest itself, equals to 200r?T,. Correspond- 
ing to it hydrostatic pressure is jH and the tangential 
component of the stress tensor on the boundary of mine 
working is 2"fH (where ^ = pg, p is rock density equal, 
in the mean, to 2.6 • l{fikg ■ m? , g is the gravity con- 
stant). From here a-r ~ \QAMPa. Assuming that shear 
hardness of rocks p, — 30GPa, we get a value for the 
shear component of the deformation tensor of the order 
of Er = cFtt/ pi = 0.00035. A square of this value has an 
order of the second invariant of deformation tensor (e.f,)^ 
that appears in free energy expansion p8)l . Taking into 
account that on the boundary of the mine working the 
first invariant of the deformation tensor is equal to zero 
(this follows from the exact problem solution in the elas- 
tic formulation), the zero condition for the determinant 
(PTjl can be written as 



/i?-4/i2(/io- 0.00035/1) =0. 



(23) 



This coupling of the theory parameters, parallel with 
attraction of the special experiment data, may be of use 
for numerical estimates. At present a complete complex 
of such experiments aimed at measurement of parameters 
/iQ, hi, /i2j M s-iid other parameters of the theory is not 
available. The consideration suggested here may serve as 
an impetus to exploratory developments in this line. 
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FIG. 4: Dependence of free energy versus average energy per microcrack during: a - a uniform tension; b - a three-axial 
unequal-component compression. 



VI. CONCLUSION 

Thus, a new mesoscopic non-equilibrium thermody- 
namic approach is developed. The approach is based 
on the thermodynamic identity associated the first and 
second \«sn of thermodynamics and having a perfect dif- 
ferential form. The introduction into the basic relation 
(HI of such internal parameters as average defect ener- 
gies and defect densities has allowed to give the phys- 
ical meaning of seeming purely mechanical problem of 
material destruction and to apply the standard methods 
closed to the equilibrium thermodynamics to the analy- 
sis of nonequilibrium processes. The identity permits to 
introduce an extended non-equilibrium state and use the 
good developed mathematical formalism of equilibrium 
and non-equilibrium thermodynamics. 

The evolution of non-equilibrium variables are de- 
scribed by a Landau-based equation set expressed 
through internal or different kinds of free energy con- 
nected by means of the Legendre transforms. The intro- 
duction of some structural viscosity-like shifts Tp^ and hu 
into the evolution equations ([5]) and ([5]) make possible 
the accordance between the different kinds of energy. 

Taking into account only one type of defects (viz., mi- 



crocracks) the description of destruction of quasi-brittle 
solids is considered by means of an expansion of free en- 
ergy. Different powers in the expansions describe thermo- 
dynamic processes of different levels. The main process 
is described by the lowest power; the processes corrected 
to it are described by higher powers of the expansions. 
Alternation of the signs is due to the Le Chatelier prin- 
ciple, namely, the processes of each next level point in 
opposite direction to the processes of the previous level. 

The consideration allows to find stable and unstable 
stationary points of free energy in dependence of sign of 
the determinant pT|) . At D > the system has the stable 
stationary state, which is separated from the unstable 
branch by energetic gap. At 13 < the system has not 
a stable stationary state and is constantly damaged with 
a creep-like scenario. A large external mechanical action 
can change the sign of the determinant from the case 
13 > to the case D < 0. Owing to this change it 
is possible to explain the origins of the high stability of 
mine workings at small depths and their instability at 
large depths. The consideration can be used for essential 
enhancement of mine working stability by the way, as 
was announced in the ref. iQl. 
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